NPS-67L177091 



LIBSASY 

TECHyxc M P??ORT 

VC. 0-1 GRADUATE 8C0&C& 
vt\ T*T i '{ , ' 'XJ^CBMA 



NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




PROGRESS REPORT ON 
AIRCRAFT FATIGUE STUDIES 

BY 

GERALD H. LINDSEY 



September 1977 



Approved for public release; distribution unlimited 
repared for : 

FEDDOCS aval Air Systems Command 

0 208.14/2:NPS-67Li77091 ode AIR-803 

ashington, DC 20361 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral I. W. Linder 
Superintendent 



Jack R. Borsting 
Provost 



The work reported herein was supported by Naval Air Systems Command, 
Washington, DC. 

Reproduction of all or part of this report is authorized. 

This report was prepared by: 



UNCLASSIFIED 



SECURITY CLASSIFICATION OF This RACE c*r >»» Omit Bnlt^td) 



REPORT DCCUUcNT AT10H PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 




1. «FO*T NUMBER 2. SOVT ACCESSION NO. 

NPS-67Li77091 


3. RECIPIENT'S CATALOG NUM3ER 


A. TITLE f»id Submit) 

Progress Report on Aircraft Fatigue Studies 


S. TYPE of report 4 PERIOD COVERED 

Interim Report 1 July 1976 ■ 
30 September 1977 




4. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR^ 

Gerald H. Lindsey 


*. CON TRACT ON GRANT HUMiERfi) 

N00019-77-WR-71012 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, CA 93940 


10. program ELEMENT. PROJECT, TASK 
AREA 4 WORK UNIT NUMBERS 

A320 3200 /136A/7R02 3-00-000 


11. CONTROLLING OFFICE NAME ANO AOORESS 

Naval Air Systems Command 
Code AIR-803 
Washington, DC 20361 


12. report date 

September 30. 1977 


11. number of paces 

31 


"JY MONITORING AGENCY name 4 AOO RESSff/ dlttcrent horn Controlling Oltice) 


IS. SECURITY CLASS, (ot thf report) 

UNCLASSIFIED 


is*, declassification / downgrading 
SCH EDULE 


i6. distribution statement (ot thta purport) 

Approved fpr public release; distribution unlimited 

* 


17. DISTRIBUTION STATEMENT (el the abctr+ct enter** in Dlock 20, It diiterent trem Report) 




16. SUPPLEWE^ 'ARY NOTES 




19. KEY WORDS (Contlrruo on rererac •!+• It ncccccany *">* f oontity by block number) 

Aircraft Fatigue 
Cumulative Damage 
Neuber T s Equation 
Residual Stress 




20. ABSTRACT (Continue on §ldo It necasoarj an* laontity by him Ji rn^b**) 

This report summarizes a year of research activity on aircraft fatigue in three 
areas: (1) The influence of minimum load levels, ground loads, order of 

loading and counting method on the damage calculation, (2) The method of cal- 
culating local stress in the plastic range at the stress concentration site from 
the recorded in-flight strain monitor and (3) The measurement of relaxation of 
residual stresses after they are produced locally at the point of stress con- 
centration. Substantial progress was made in all three areas and especially 




DD , j°" M 73 1473 COITION OF 1 NOV *» IS OBSOLETE 

(Page 1) S/N 0 102-014- 4601 | ^ SECURITY CL AS JI F 1 C ATION OF THIS FAOt (*»•" 3ml. 



s C CURITY CLASSIFICATION OF THIS RAO* C»>*" 0»*» 



ABSTRACT 



This report summarizes a year of research activity on aircraft fatigue in 
three areas: (1) The influence of minimum load levels, ground loads, order 

of loading and counting method on the damage calculation, (2) The method of 
calculating local stress in the plastic range at the stress concentration 
site from the recorded in-flight strain monitor and (3) The measurement of 
relaxation of residual stresses after they are produced locally at the point 
of stress concentration. Substantial progress was made in all three areas 
and especially the first two. 



iii 



TABLE OF CONTENTS 



Page No . 

I. INTRODUCTION AND BACKGROUND 1 

II. DAMAGE CALCULATION 4 

III. LOCAL STRESS CALCULATION 12 

IV. RESIDUAL STRESS RELAXATION 17 

V. SUMMARY AND CONCLUSIONS 21 

VI. REFERENCES 29 

VII. DISTRIBUTION LISTS 31 



v 



I. INTRODUCTION AND BACKGROUND 



Early Developments 

Commencing three years ago, a feasibility study was undertaken to determine 
if an effective aircraft fatigue monitoring instrument could be conceived on 
the basis of microprocessors* The goal was to design a small, lightweight 
package, which could record maximum and minimum loads sequentially in a form 
conducive to automated data processing. The first unit was designed and 
fabricated by LT David Vidrine (1) and LT James Sturgis (2). This first pro- 
totype was physically larger than would have resulted for an optimum design 
solely for fatigue monitoring in order to give it testing, check-out, and 
development capability. It was built around the Intel 8008 CPU chip, and had 
16 channel capacity via a multiplexer. It digitized the incoming load infor- 
mation from strain gages, accelerometers, or other transducers and then 
searched these signals for maxima and mimima. Upon finding them, the digital 
value of each extremum point was stored in memory in chronological order until 
the memory was filled, at which point, the microprocessor caused the tape 
recorder to be activated, and the memory to be transferred to the tape cassette. 

This unit was laboratory tested with sinusoidal signals from signal 
generators with results as shown in Fig. 1. After success was obtained with 
known signals, the unit was flight tested in an NPS leased Cessna 310. Typi- 
cal results are shown in Figure 2. 

Second Prototype 

In those budding days of microprocessors, the state of the art was 
advancing rapidly as it is today, and by the time the first prototype was 
built, new developments had been perfected. LT Wesley Stanfield (3) began 
work on a second prototype that was directed toward flight test in a high 
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performance aircraft. It was based on the 8008 CPU, but incorporated ad- 
vancements in programming capability and reduction in the physical size of the 
components. It was designed with 8 channel capacity and had the physical 
dimensions shown in Figure 3. The final version of the software was generated 
by LT Charles Butler (4) , who also generated the programs for data reduction 
from the instrument through a Hewlitt Packard 9830, desk- top computer. 

This second prototye was successfully laboratory tested and demonstrated 

for Dr. David Weiss of the Naval Air Development Center (NADC) . The instrument 
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was able to take data from signal generators up to 5H , pick out maxima and 
minima and record them on tape. Using the HP 9830, the voltage of the maxima 
and minima are listed in tabular form in chronological sequence. 

Final Results 

This concluded the feasibility study at NPS, and the Naval Air Systems 
Command, through Bob Weinberger, adopted the concept of an in-flight micro- 
processor instrument to monitor fatigue as well as other flight and engine 
parameters. The F-18 was designated as the first Navy aircraft to incorporate 
the concept, and all future aircraft will use state-of-the-art versions of this 
fatigue monitoring concept. 

The second prototype instrument was then to be transitioned to Dr. Weiss 
at NADC for further experimentation while new prototypes were designed and 
packaged to meet military specifications. A malfunction in the analog to 
digital conversion unit was experienced at this time. Steps taken to correct 
it finally required a new unit from an alternate manufacturer, which in turn 
precipitated changes in the wiring and associated changes in the software. In 
the meantime, further developments in microprocessor technology have made 
significant investments of time in the second prototype unit unwarranted. 
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Future Possibilities 



In the original prototyping, a tape recorder was used as the final 
reservoir of the data. This was a temporary arrangement because of the limited 
conditions under which successful tape recorder operations can occur. Large 
temperature and humidity variations can cause the capstans to freeze, and 
recording in a high g environment may not be reliable at all; however, there 
is an alternative emerging. 

Two specific possibilities have recently come forth: one is a CMOS storage 
cell, which requires so little power that small batteries can sustain them 
into the time ranges of weeks that would be needed for actual operations. 

The second is Texas Instrument’s bubble memory, which at this time is not 
completely reliable but which has promise to fill the needs of the system in 
the not too distant future. by providing large quantities of nonvolatile memory 
capacity. 
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II. DAMAGE CALCULATION 



In light of the Navy commitment to employ microprocessor units as fatigue 
monitors, it now becomes appropriate to look forward to the type of data that 
will be available from them and how it can be used to calculate fatigue 
damage. 

Two General Approaches 

The Navy for many years has used a cumulative damage approach to fatigue, 
based upon Miner T s Law or some variation of it, which has now been developed 
into a rather complex calculation that considers various methods of counting 
the cycles, plastic flow around the point of stress concentration, residual 
stresses and their relaxation, etc. On the other hand the Air Force is pur- 
suing fatigue calculations based upon crack propagation. An initial flaw of 
.050 in. is assummed to exist at each critical stress riser location, and the 
number of cycles required to grow it to a critical size constitutes the 
fatigue life. 

Two large computer programs have been imported and made operational on 
the NPS IBM 360. Potter’s (5) cumulative damage program uses Miner’s Law and 
is presently operational. Engels (6) program, "Cracks II, " employ the crack 
propagation method and is now being checked out with example problems provided 
by the author. Check runs are planned in the future to give comparisons of 
fatigue lives for the same conditions and to investigate the affect of cycle 
minimums in the load spectrum. The calculations will be for an actual air- 
plane fatigue critical point on the A-7E, 

Potter Program 

LT Scott Atkinson (7) has investigated the influence of counting method 
and type of block loading on damage. The load spectrum used in the analysis 
was the spectrum which is reproduced in Table I, from MIL SPEC A8866, and the 
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TABLE I 



FREQUENCY OF MANEUVER LOADS 

NUMBER OF TIMES PER THOUSAND HOURS THAT 
LOAD FACTOR IS EXPERIENCED 



PERCENT OF MAXIMUM (POSITIVE) FLIGHT MANEUVER LOAD 

SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 



35 


17000 


45 


9500 


55 


6500 


65 


4500 


75 


2500 


85 


1360 


95 


440 


105 


150 


115 


40 


125 


16 



TOTAL 42006 



PERCENT OF MAXIMUM (NEGATIVE) FLIGHT MANEUVER LOAD 

SYMMETRICAL LIMIT LOAD FACTOR SPECTRUM A 



0 


500 


10 


200 


20 


100 


30 


60 


40 


35 


50 


30 


60 


25 


70 


20 


80 


15 


90 


10 


100 


5 


110 


3 



TOTAL 1003 
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loads are assumed to be applied at wing station 32 on an A-7E wing (8). Load 
data were input as block loads and also as single randomized cycles. The 
Miner's Law summation of fatigue damage for this loading spectrum is .1300 
per 1000 flight hours. Results obtained using this computer program are 
given in Tables II, and III. It is important to notice that when fatigue 
damage is calculated without using the range-pair counting method, it makes 
only a minor difference on the total fatigue damage when the sequence of 
loading is changed or when negative loads are included. 

Block loads were arranged in L0-HI, HI-LO, and HI-LO-HI sequences. The 
results indicate the LO-HI sequence to be more damaging than the HI-LO 
sequence, with damage due to the HI-LO-HI sequence falling between the other 
two block load sequences. In the HI-LO sequence the HI loads leave a local 
residual compressive stress in the material, and the LO loads do less damage. 

In the case of randomly generated load sequences, identical starter 
integers were used in order to evaluate the effect on the damage of inclusion 
of MIL spectrum A negative loads and ground cycles. The more negative loads 
that are included, the worse is the fatigue damage. The increase in damage 
with more negative loads is caused by the negative loads decreasing the local 
residual compressive stress caused by positive loads into the plastic zone. 

This permits the local positive elastic stresses to be more damaging. 

Minimum Load Levels 

Including the negative values of stress resulting from landing loads, 
which through the range-pair counting method is^ tantamount to appropriately 
accounting for the ground-air-ground cycle, makes considerable difference in 
damage accumulated even for this high performance airplane. This is, of course 
contrary to many statements in the literature that GAG cycles are only of 
significance in transport aircraft, and this conclusion is proper when 
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TABLE II 



FATIGUE DAMAGE DUE TO POSITIVE MILSPEC A 
BLOCK LOADING, NORMALIZED TO 1000 HOURS 



SEQUENCE OF 
BLOCK LOADS 


11% LIMIT LOAD 
USING RPCM* 


MIN. 


11% LIMIT LOAD ; 
NOT USING RFCM* 


LO-HI 


.16780102 




.16780108 


HI-LO 


.069243014 




.069243014 


HI-LO— HI 


.10864216 

* RANGE-PAIR COUNTING METHOD 
TABLE III 

FATIGUE DAMAGE DUE TO MILSPEC A 


.10839599 


RANDOMIZED 


LOAD INPUTS, USING RPCM, 


NORMALIZED TO 


1000 HOURS 


INTEGERS USED 


TO START RANDU 


MINIMUM 


LOAD VALUES 






11 PERCENT 


NEGATIVE 


GROUND 




LIMIT LOAD 


MILSPEC A 


CYCLES 


83745,54711,54487 


.039226338 


.046267733 


.060848035 


13547,66549,7 


.095124505 


.085063167 


.10429338 


9,583,4777 


.043321513 


.043640956 


.036781987 


48621,3,491 


.018672124 


.034077277 


.035881512 


73,559,1001 


.031819174 


.042006299 


.037130679 


357,833,1 


.037665060 


.049336702 


.050973855 
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residual stresses are not being accounted for. However the negative loads 
cause a reduction in the protective residual stress at the stress concentration 
site and an associated increase in damage as shown by the third column of 
Table III. 

The entire question of minimum loads is raised at this point. Since 
minimum loads are not recorded at each cycle on the accelerometer counters, 
it has always been assumed that the minimum of each cycle was lg. , which is 
intuitively thought to be unrealistic. However from the data in Table III, 
the minimums are important and that information will be recorded and available 
for use from the new microprocessor fatigue monitoring instruments. In the 
meantime, flight testing on current aircraft with the same instruments can 
be used to generate a spectrum of minimum load levels for each cycle. This 
spectrum can then be randomized and used in the program in the same way that 
the 11% loads and the negative MILSPEC A loads were used in this study. 

Random Load Histories 

Figure 4 illustrates the pattern of the individual load randomization 
process and portrays the first two flight hours of the aircraft. The range of 
values in the columns of Table III are caused by such variations in the load 
randomization pattern. 

The fatigue damage resulting from randomization of loads on a cycle basis, 
as shown in Table III, is significantly different from the damage due to the 
HI-LO-HI block loads illustrated in Table II, which has been accepted by many 
as being representative of random loading. Eor instance, Vought Aeronautics 
Division recently conducted a fatigue study using randomized flights where 
each flight was based on a mission profile specified by the Air Force. They 
found fatigue damage resulting from the randomized flights was in close agree- 
ment with fatigue damage resulting from HI-LO-Hi block loading. 
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Table IV has been compiled to show the influence of not using the range- 
pair counting method, van Dijk ! s (9) arguments are so compelling concerning 
the advisability of using the range-pair method and these results show that it 
should be used although it is considerably more expensive to run. The range- 
pair counting method relates the load sequence to the stable, cyclic, stress- 
strain behavior of a material by coupling up strain ranges so that each 
combination of half cycles will produce a complete, closed stress-strain hyste- 
resis loop in order to be counted as a cycle. 

The range-pair counting method is illustrated in Figure 5. It counts a 
strain range as a cycle if it can be paired with a subsequent straining of 
equal magnitude in the opposite direction. For a complicated load history, 
some of the ranges counted as cycles will be simple ranges, such as 2-3, 
during which the strain does not change direction, but others, such as 1-8, 
will be interrupted by smaller ranges, which will also be counted as cycles. 

In Figure 5, ranges are marked with solid lines and the paired ranges with 
dashed lines. 

Each peak is taken in order as the initial peak of a range, except that 
a peak is skipped if the part of the history immediately following it has 
already been paired with a previously counted range. If the initial peak of 
a range is a minimum, a cycle is counted between this minimum and the most 
positive maximum which occurs before the strain becomes more negative than 
the initial peak of the range. For example, in Figure 5 a cycle is counted 
between peak 1 and peak 8, peak 8 being the most positive maximum that occurs 
before the strain becomes more negative than peak 1. If the initial peak of 
a range is a maximum, a cycle is counted between this maximum and the most 
negative minimum* which occurs before the strain oecomes more positive than 
the initial peak of the range. For example, in Figure 5 a cycle is counted 
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between peak 2 and peak 3, peak 3 being the most negative minimum before the 
strain becomes more positive than peak 2. Each range that is counted is 
paired with the next straining of equal magnitude in the opposite direction, 
explaining why complete cycle rather than half cycle counts are made. For 
example, in Figure 5 part of the range between peaks 8 and 9 is paired with 
the range counted between peaks 1 and 8. 
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TABLE IV 



FATIGUE DAMAGE DUE TO MILSPEC A RANDOMIZED LOAD 



INPUTS , 


NOT USING RPCM, 


NORMALIZED TO 1000 


HOURS 


INTEGERS USED 


TO START RANDU 


MINIMUM LOAD VALUES 






11 PERCENT 


NEGATIVE 


GROUND 




LIMIT LOAD 


MILSPEC A 


CYCLES 


83745,54711,54487 


.040602274 


.047398917 


.049124956 


13547,66549,7 


.040137805 


.047727041 


.049612150 


9,583,4777 


.041204430 


.048346408 


.049540661 


48621,3,491 


.041160211 


.049230121 


.050957575 


73,559,1001 


.040792227 


.047043338 


.049311966 


357,833,1 


.041056350 


.049220584 


.050481893 
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III. LOCAL STRESS CALCULATION 



There are two aspects of this program that require further investigation, 
and the first of these, local stress calculation will be discussed in this 
section. By local stress is meant the stress at the site of the stress 
concentration. The transfer function between the local stress and the far-field 
strain of a skin in uniaxial tension must be a function of the geometry of the 
notch or hole, the local strain, and the type of material. 



a = f (E,p , a, e, s, e) 



( 1 ) 



where 



E = Young’s Modulus 

p = Geometric parameter - radius of curvature 
of the notch or hole. 

a = Geometric parameter - size of the notch or 
hole 

e = Local strain 
s = Far-field stress 
e = Far-field strain 



To nondimensionalize this set of variables, stress concentration factors 
are used. For instance any notch can be fit with an ellipse as proposed by 
Inglis (11), and the stress concentration factor for an ellipse in the elastic 
range is 

K t i f = [1 + 2 /f] (2) 

Defining dimensionless ratios for local stress and strain in the plastic 
range in the natural way. 
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the transfer function can be simplified to 



= / CK c> 19 



( 4 ) 



Neuber (12) proposed a function to describe this relationship in the form 
of a geometric mean of the stress and strain concentration factors, 

2 



K K = K_ 
as T 



(5) 



This was evaluated on 7075-T651 experimentally by Garske (13) by using large 
notched plates so that strain gages could be placed at the stress concentration 
for reliable local strain readings. Four plates were studied with = 3.92 , 
3.33 , 2.31 and 1.63 . Neuber 1 s relationship was found to be in error in 
predicting local stresses and strains by as much as 24% in the higher plastic 
range in all four sheets. 

From the appearance of the data in graphical plots and guided by Neuber f s 
original proposal, a slightly generalized form of the power law was proposed. 

K X K 7 = A K Z (6) 

a e T 



To properly maintain the elastic limit condition where = K = K , A = 1 
and z = x + y . 





(7) 



Without loss of generality, 
Rearranging, 



y can be set to one for later ease in calculation. 




( 8 ) 
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A curve fit of the stress and strain concentration in the plastic range was 
made of the data for plate //3 for which K = 2.51. For a correlation 
coefficient of .90, the result was 




(9) 



The fact that this came very close to the proper elastic limit constant of 
unity was very encouraging that the general form is correct. Whether there is 
something fundamental in the form of the equation and 1.24 should be 5/4 is not 
known. 

It has the general appearance of an energy density, but after several 
attempts to derive it from fundamental principles failed, it was taken as an 
empirical law as it is. However to keep the appropriate and consistent elastic 
form, the equation should be 




( 10 ) 



In further calculations this equation (10) will be used. 

To carry out the calculations, the stress-strain relationship was formulate 
from experimental data on 7075-T651 in the classical Ramburg-Osgood form by 
a power law curve fit of the plastic stress-strain data points. 



£ 









a 


-L. 


0 


,10.22xl0 6 . 


T 


99815^ 



22.678 



(ID 



This expression fit the nonlinear portion of the uniaxial data with a correla- 
tion coefficient of .98. 

The final expression of local stress as a function of far-field strain can 
now be formulated from the rearranged geometric equation (10) . 
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( 12 ) 



1.244 



(f) ’ K ' 



2.244 



Since s = Ee , 



(t) 



1.244 



= O^e) 



2.244 



(13) 



Inserting the stress-strain expression for e , 



i 2. 244 



,23.922 



VlO. 22x10' 



VI. 27x10' 



- ( V ) 2 - 244 - 0 



(14) 



Solutions of this equation for stress over a range of values of far-field 
strains gave maximum variations of 4% from the test data for the three 
different K^'s tested. 

With the local stress determined accurately, the residual stress can be 
determined accurately by assuming elastic unloading. For instance, if the 
fatigue monitor returned with a sequence of recorded strains as shown in 
Table 1, the corresponding local stresses could be calculated by the techniques 
that have been developed here. 

TABLE I 
K t = 1.63 

Cycle Far-field Far-field Local Stress 





strain 


stress (psi) 


(psi) 


1st Peak 


.005308 


54,248 


77,163 


1st Valley 


.000000 


0 


- 11,260 


2nd Peak 


.001455 


14,870 


12,960 


2nd Valley 


.000000 


0 


- 11,260 
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The first peak local stress comes from equation (14) using the given far-field 
strain and . The far-field stress producing it is s = E (*005308) = 54,24; 
The first valley stress is calculated by removing the far-field stress of 
54,248 psi, which corresponds to locally removing Ag = K^S = (1. 63) (54 ,248) = 
88,424 psi. Subtracting 88,424 psi from the peak stress of 77,163 psi gives 
the residual stress of - 11,260 psi after the first cycle. The second cycle 
strain of .001455 produces a stress increase from equation (14) of 24,226 psi. 
However since there was a residual of - 11,260 psi, the second cycle peak 
stress will be 12,966 psi. Upon unloading, since there has been no plastic 
yielding, the structure returns to the same residual stress state as before. 

The first cycle yielding has been of considerable benefit in protecting the 
structure during the second cycle with a compressive residual stress, which 
reduces the fatigue damage accumulated. This can be repeated throughout a 
randomized or an actually recorded load history where the range-pair counting 
method has been incorporated. 
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IV. RESIDUAL STRESS RELAXATION 



The second of the two aspects of the program that requires investigation is 
the relaxation of the residual stress created by loads into the plastic. In the 
true sense of the word, stress relaxation is stress decay under constant strain, 
but for the present the term stress relaxation will be used more loosely to mean 
stress reduction with time and will be descriptive of residual stress behavior 
with time. 

When structures containing stress concentrations are loaded sufficiently to 
produce local plastic straining in the region of high stress, residual compres- 
sive stresses remain in that region upon unloading of the structure, as has been 
described. The residual stress and strain are highly localized disturbances in 
the elastic continuum, and since they are disturbances, they are potentially 
unstable. Given sufficient conditions, they will tend to be reduced and will 
relax toward lower values during succeeding cycles until the next overload occurs. 
Potter (10) has proposed that the residual stress be decomposed into an equilib- 
rium portion and a transient portion and suggested an exponential to model the 
transient behavior, thinking of it as being similar to a transient in a critically 
damped dynamic system. The overload triggers the impulse and controls its height, 
and following the impulse, there is an asymptotic return to an equilibrium state. 
Potter* s model for residual stress, o R , could be expressed as follows: 

a R = CT REQ + a RTRANS (15) 



In terms of the proposed exponential, 



ct_ = a + Ae 
R REQ 



bN 



( 16 ) 



N = Number of cycles 

where it is assumed that the transient residual stress decreases with cycles 
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A and b are evaluated from the initial conditions and the end conditions. Xf 
a R = a RQ at N = 0, then 



A = 



a R0 ~ a REQ 



(17) 



If N= N__ is the point when the transient residual stress is 90% relaxed; i.e. 
EQ 

= a„ - .9 (a__ - c__,_) 

R R RO REQ 



b 



In 0.1 




(18) 



this gives for 



R’ 



a R a REQ + ( °R0 



* ,ln 0.1 N . 

°R£Q )e( H eq ) 



(19) 



or 



= a. 



REQ 



+ (a 



RO 



a REQ )e 



- 2.3 N 



N. 



EQ 



( 20 ) 



Experimental tests on uniaxial specimens were run by Bentley (14) to eval- 
uate the undetermined parameters in this model and study this general behavior. 
The first area of interest was in the mechanism of the substructural recovery, b? 
relaxation of the lattice structure, that accounts for the stress decay toward 
a stabilized stress level. Cyclic relaxation tests on 7075-T7651 into the plast: 
zone between fixed strain levels revealed no detectable stress relaxation. Con- 
sistently every relaxation test showed the reloading stress— strain curve rejoin* 
at the point of departure from the original stress— strain curve for the cyclic 
unloading, which is evidence of zero relaxation. However, another recovery 
mechanism, anelastic behavior, was measured, and it was concluded that it was 
the primary element producing residual stress decay. 

Anelastic behavior is illustrated in Figure 6, where it is seen that all of 
the plastic strain present at unloading is not permanent set, but some of it is 
recovered. After the load was removed, the e continued to decrease exponen- 
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tially with time for periods of the order of one hour. Using the same form of 
the model as we did for stress. 



and 



e = e 



equil. 



+ £ 



'transient 



( 21 ) 



s - (e - e )exp(-2. 303 t/t ) (22) 

trans p eq F eq 

At a point of stress concentration in the structure, the material under 
highest stress may be thought of as a miniature uniaxial specimen as shown in 
Figure 7. After plastic straining at this location, the specimen is longer 
than it was originally, and the surrounding elastic material is trying to com- 
press it back to its original length. This puts the element in compression and 
is where the residual stress arises from. However, as the plastic strain 
decreases through anelastic recovery the mismatch is gradually reduced and not 
as great as at first. The surrounding material, consequently, does not exert 
as great a compressive stress as it did originally. The transient portion is 
given by 



<7_ = ( a r»^” a r* )exp(-2.303 t/t ) 

Rtrans R0 Req r eq 



(23) 



This process was found by Bentley to be time dependent rather than cycle 
dependent. The anelastic recovery was measured on one specimen undergoing a 
single load cycle, which was then left to sit. It was then measured on a specimen 
which was cycled continuously. The anelastic recovery was the same in each case. 
The residual stress decay is a function of time and not cycles. 

The total expression for residual stress is 



a D = a + (a a ) exp (-2. 303 t/t ) 
R Req RO - Req eq 



(24) 



If a is expressed as a percentage of o , the anelastic recovery tests 
Keq RO 

can be used to evaluate the unknown parameters in the model. The a is pro- 

RO 
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duced by the mismatch in the strains at the edge of the notch 



RO 



= k e 



( 25 ) 



The equilibrium stress, a , is produced by the mismatch in strains remaining 
after time. 



Req 



= k e 



eq 



= a 



RO 



e 

P 



(26) 



Using these in the expression for yields 



a R CT R0 



, E 

_e£ 

£ 

P 



) + i 1 - T 1 ) exp (■ 2 ' 303 c/t eq) 



(27) 



Tests were run at two different strain levels in the plastic zone with results 
as given in the following table: 

TABLE II 



£ 


£ 


e /e 


t 


P 


eq 


eq p • 


eq 


.00540 


.005033 


.932 


53.49 min 


.005046 


.004784 


.948 


50.05 min 



A broader range of strain levels in the plastic zone would need to be run to deter 

mine if e / e is a constant and if t is a material constant, 
eq p eq 

If this were the case, the residual stress expression would become 

a D = [0.94 - 0.06 ex P 0.0919O] ' (28) 

K RO 

where a D is found as described in the last section. 

Certainly these few tests cannot be relied upon to accept this expression 
as factual for 7075-T651; however two important points have been forthcoming. 
First, the stress decay was dependent upon time and not upon cycles, and secondly 
it was small in magnitude. All of this will be probed more thoroughly and deeply 
during the next year . 
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V. SUMMARY AND CONCLUSIONS 



A fatigue monitoring instrument based upon a microcomputer has been demon- 
strated in laboratory and preliminary flight test to be feasible and very desir- 
able* Due to its adoption for future aircraft by NAMAIR, investigations have been 
begun for effective use of the far-field strain data recorded by it to compute 
fatigue damage. 

Three major areas have been undertaken in this study: (1) A sophisticated 

computer program has been used to investigate the effects of the minimum portion 
of each load cycle on the accumulated fatigue damage. (2) Analytical and expe- 
rimental investigations have been made of the use of Neuber's equation in calcu- 
lating local stress from far-field strain with the result that a better equation 
was developed. (3) An experimental investigation of the relaxation of residual 
stresses produced at sites of stress concentration was undertaken with some 
success but with other work to be done before these results are conclusive. 
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Figure 1 Laboratory Test Results of the 
Fatigue Monitor 
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Figure 3 



MIDAS Aircraft Fatigue Monitor 
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Strain 



Figure 5 




Illustration of Range-Pair Counting Method 
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Figure 6 An Illustration of Anelastic Behavior 
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Figure 7 Uniaxial Representation of Residual Stress 
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